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Under the assignment of D s (3040) as a radially excited p-wave cs state with J p = 1 + , we compute 
the P s —¥ D s (3040) form factors within the covariant light-front quark model. Two classification 
schemes for the p-wave cs meson are adopted. We also use our results to predict the branching 
ratios (BRs), polarization fractions, and angular asymmetries in semileptonic P B — > D s (3040)te(^ = 
e,fi, t). The BRs are found at the order of 10~ 3 for £ = e,fi, and 10~ J for I — r. We find that 
the polarization fractions and angular asymmetries could be useful to pin down the ambiguities in 
theoretical prescriptions for D s (3040). In addition, we investigate the nonleptonic P s — > D s (3040)71'/ 
decays under the factorization method, where M denotes a charged pseudoscalar or a vector meson. 
The BRs of P 3 -> D s (3040)p and P s -> D S (3040)D* reach the order of 10" 3 , while the other 
channels are typically smaller by 1-2 orders. 

PACS numbers: 12.39.Ki; 13.20.He; 13.25.Hw 



I. INTRODUCTION 



The observations of a number of new resonances, such as the narrow states D S (2317) in the D+ir° final state and 
D s (2460) in the D*ir° and D s j channel [HQ, have brought great opportunities for our understanding of nonpertur- 
bative QCD in the charm mass region. They also initiate tremendous efforts on the exploration of exotic hadrons 
beyond the naive quark model and revive great interests in the study of the cs spectrum Q . 

In the heavy quark limit, the heavy quark will decouple from the light degrees of freedom and act as a static color 
source. Strong interactions will be independent of the heavy flavor and spin. In this case, heavy-light mesons, the 
cigenstates of the QCD Lagrangian in the heavy quark limit, can be labeled according to the total angular momentum 
s; of the light degrees of freedom. Heavy mesons with the same s; but different spin orientations of the heavy quark 
fill into one doublet. Accordingly, heavy mesons can be classified by these doublets, characterized by s; instead of the 
commonly-used quantum numbers, 2S+1 Lj. 

The classification scheme based on the heavy quark symmetry applies to the D s j mesons. For instance, £> s (1969) 
and Z?*(2112) can be arranged into the lowest-lying doublet with Jf = (O^,! - )^- Among the p-wave cs states, 
D s0 (2317) and Z? s i(2460) belong to the doublet with J? = (0+,l + )i /2 , and doublet j£ = (l+,2+) 3/2 can be filled 
by -D s i(2536) and £> S 2(2573). Despite some controversies over their internal structures, these mesons have been well 
established in experiment. 

Several other resonant signals have also been observed recently on the B factories and other facilities. The SELEX 
collaboration reported one narrow state D s ,/(2632) in 2004 3, which however, has not been confirmed so far by other 
experiments. D s (2710) is announced in B meson decays by the Belle collaboration [f| and also confirmed by the 
BaBar results jg|. The D s (2860) is discovered in the DK final state by BaBar Recently, the BaBar collaboration 
has reported one new D s j state in the study of the e + e~ annihilations 0. The mass and decay width of this broad 
state are 

m = (3044 ± 8jJ° )MeV, T = (239 ± 35±^)McV. (1) 

Its quantum numbers could be J p = 0~, 1 + , 2~..., taking into account its signal in D*K instead of DK. The angular 
analysis is not available because of the limited statistics. 

During the past few years, some progresses have been made in the classification of these states and understanding 
of their peculiar properties in the production or decay processes. Towards this direction, if the 3 _ assignment of 
£> s (2710) @-[l3] is confirmed by the data (see Ref. [U for an alternative scenario), the D s (3040) would be unfavored 
to be 2 _ since the mass inversion seems unlikely. On the other hand, the 1 + assignment is favored in several aspects. 
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Firstly, this explanation is supported by the study of the spectroscopy of excited D s mesons I12h14| , though predictions 

m Refs. dEEai 

are lower than the experimental value by roughly 200 MeV. Reference [10( investigates these new 
D s j mesons in a semi-classic flux tube model and concludes that D s (3040) is compatible with the 1 + assignment. 



The authors of Ref. I7[ have studied the strong decays of D s (3040) with four different assignments, and find that the 



decay width with J = 1 + is also consistent with the value given by BaBar. The strong decays of Z? s (3040) are also 
investigated in the 3 Po constituent quark model [l8| and chiral quark model [l9j . respectively. 

The verification of the above new states can be via not only the analysis of their decay properties, which have been 
widely discussed in the literature, but also via their productions, which so far have not been investigated in detail. 
The production processes of the above mentioned new states in B s decays would be an ideal way of probing their 
properties. Productions of D s (2317) and D s (2460) in the B s decays have received some theoretical attentions |22T[26| . 
On the experimental side, a large number of B s events will be collected by the LHCb [2(| and forthcoming Super 
B factory experiment [2l[. From the theoretical viewpoint, the degrees of freedom over the nib scale in B s decays 
can be computed in the perturbation theory and the evolution between the mw and nib can be organized using the 
renormalization group. The low-energy effects, which are parameterized into hadronic form factors, will then probe 
the structures of the D s j mesons. 

To proceed, we will study the B s — > D s j form factors and the production rates of D s j mesons in semileptonic and 
nonlcptonic B s decays in the framework of a covariant light- front quark model (LFQM) [27| . This approach is suitable 
for handling the heavy- light mesons and has been successfully applied to various processes p8l-[33j. In particular, we 
will focus on the £> s (3040) with its assignment as a 2P state. The B s — > D* transition form factors will be extracted 
as a byproduct. 

This paper is organized as follows. In the next section, an introduction to the light-front quark model and expressions 
for the transition form factors between the B s and (D*, D S (3QAQ)) mesons are presented. In Sec. Ill, we compute 
the form factors under two different assignments for the D s (3040). With the form factors at hand, we will explore in 
Sec. IV the semileptonic B s decays and the nonlcptonic B s — > D s (3040) decays under the factorization assumption. 
The last section contains our conclusions. 



II. FORM FACTORS IN THE COVARIANT LFQM 

In the effective electroweak Hamiltonian responsible for B s — > D s jlv [34| 



G 



H cS = ^ c6 [c Tm (1-75)&][W-75H (2) 

Gf and V c b are the Fermi constant and the Cabibbo-Kobayashi-Maskawa (CKM) matrix element, respectively. The 
leptonic part can be directly calculated using perturbation theory. The residual part contains hadronic effects and 
will be incorporated into the transition form factors as 



(D* s (P",e"*)\V^\B s (P')) = ^ ^ a pe"*»P a q eV B ° D °(q 2 ), 

rriB 3 + m D , 

(D* 8 (P",e"*)\A fi \B s (P>)) = 8 {K+m D .)^f c; fe 2 )- £ "\ '- P P,A^ D ° (q 2 ) 
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q ,[A^ D H q 2 )-A^ D Hq 2 )]}, 



(D sl (P",e"*)\A fl \B s (P')) = e liya pe"*»P a qPA B ° D ^(q 2 ), 



(D sl (P",e"*)K\B s (P>)) = -{{m Bs m^J^F^ ( g 2 ) - -1 — — - P,V 2 B ^(q 2 ) 

- 2mDsl £ -^q,[V 3 B ^(q 2 ) - V B ^(q 2 )]}, (3) 



where P = P' + P", q = P' — P" , and the convention £0123 = 1 is adopted. The vector and axial-vector currents are 
defined as c^f^b and 07^756. We have adopted the 1 + assignment, denoted as D s i hereafter, for the f s (3040) meson. 
In the subsequent analysis two kinds of axial- vectors { 2S+l Lj = 3 Pi or will be considered, and for short we will 
also abbreviate them as 3 A and 1 A, respectively. To smear the singularity at q 2 = 0, we obtain the constraints on the 
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FIG. 1: Feynman diagram for a transition form factor, where the cross symbol denotes the electroweak vertex. 



form factors Af sD: (0) = Ag sD ' s (0) and V^^O) = V BsDsl (0), where 

vf-°»tf) = - - + v«-°- (,'). (5) 

2m Dsl 2m Dsl 

To be specific, we will compute these transition form factors by considering the general P — > V and P — > A 
transitions, where A denotes either 3 Pi or 1 P\ state. The LFQM is based on the light- front field theory, where 
the plus component of the gluon degrees of freedom vanishes in the light-cone gauge. In this framework, it is 
convenient to use the light- front decomposition of the momentum P' = (P'~ , P' + , Pj_), with P /=t = P'° ± P' 3 , so that 
p/2 _ p<+p>- _ p<2 The incoming (outgoing) meson has the momentum of P' = p[ + P2 (P" = p'{ + P2) and mass 

of M' (M"). The quark and antiquark inside the incoming (outgoing) meson have the mass m'/ and 7712. Their 
momenta are denoted as p'S and p-2, respectively. In particular, these momenta can be expressed in terms of the 
internal variables (xi,p'j_) 

Via = x i,* p ' + > P'i,2± = xi,2P±±p'±, (6) 
with x± + X2 = 1. Using these internal variables, one can define some useful quantities for the incoming meson 



where can be interpreted as the energy of the quark or the antiquark, and M is the kinetic invariant mass of the 
meson system. One advantage of these internal quantities is that they arc Lorentz invariant. The definition of the 
internal quantities for the outgoing meson is similar. To formulate the amplitude for the transition form factor, we 
also require the Feynman rules for the meson-quark-antiquark vertices («r' fl/ ) 

iT'p = H' pl5 , 

iV' v = iH'vYi^ - pprVi -P2)^\, 

iT' 3A = iHLjJr/p+up-ty'i-pa)^, 
yvs A 

iT[ A = i Hi A [-^(p' 1 -p 2 )^ 5 . (8) 

In the case of the outgoing meson, we shall use 2(701^70) for the relevant vertices. 

The lowest order contribution to a form factor is depicted in Fig. [TJ in which the cross label in the diagram denotes 
the V, A transition vertex. In B s to D s j decays, Pi(p'{) is the momentum of the bottom (charm) quark, while P2 
is the momentum of the antiquark. Physical quantities, e.g. decay constants and form factors, can be expressed in 
terms of Feynman momentum loop integrals which are manifestly covariant. As an example, we consider the P — > V 
transition 

E PV _ -3 N c f ,4 / H'pjlHy) oPV (q s 

B » - 1 (2^)4 J d ^ N [N'{N 2 ^ ' ( ' 
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m'i"^ 2 and N 2 — p\ — m| arise from the quark propagators and N c = 3 is the color factor. The 
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function S^f is from the Lorentz contraction of the propagators and the vertex functions 
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(10) 



The form factors will be obtained by evaluating the above formulas. In practice, we perform the integration over 
the minus component within the contour method. If the covariant vertex functions are not singular when perform the 
integration, the amplitude will pick up the singularities in the antiquark propagator. This manipulation then leads to 
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with p" = p'j^ — x 2 q± ■ The explicit forms of h' M and w' M used in this work are given by (27l . [25 
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(13) 



where tp' and ^ are the light-front wave functions for the s-wave and p-pave mesons, respectively. 

One difference between the conventional LFQM [35|-f38| and the covariant LFQM lies in the treatment of the 
constituent quarks. In the conventional LFQM, the constituent quarks are required to be on mass shell. Physical 
quantities can be extracted from the plus component of the current matrix elements. However, this framework suffers 
from the non-covariance problem arising from the missing zero-mode contributions. This drawback is resolved by 
including the so-called Z-diagram. In the covariant LFQM the quarks and antiquarks are off-shell and the physical 
quantities are written as the integration over the 4-momentum of the antiquark. The conventional LFQM will be 
recovered after the integration over the minus component, where the antiquark is set on-shell. The Lorentz covariancc 
would be lost again as it receives additional contributions proportional to the light-like four vector a) = (0,2,0j_). 
This is spurious since the vector never appears in the definition of the form factors. The advantage of the covariant 
LFQM is that we can directly handle the non-covariant terms instead of compute the Z-diagram. Specifically, the 
spurious terms can be eliminated by performing the p~ integration in a proper way and the corresponding rules are 
derived in Refs. [13, H^j. The above manipulation results in the expression for P —> V form factors: 
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The functions Z 2 , A^\ A^\ and A^ arc given in the Appendix. 

The computation of B s — > Z? s (3040) transition is analogous, and in fact, the B — s- 3 A form factors are related to 
the B — > V ones. In particular, the auxiliary form factors are given as 



1 A (q 2 ) = f(q 2 ) with (m'/ -777'/, h\ 



l3 A , w'v 



q A {q 2 ) = g(q 2 ) with (mi -> -m'/, -> ^' A , w£ -> utf A ), 

cl A (g 2 ) = a±( g 2 ) with (777'/ -> -m'/, -> ft.'4, -> 



(16) 



where the replacement of m'/ — > —777'/ does not apply to m'/ in w" and ft.", as it arises from the propagators and 
quark-antiquark-meson coupling vertex 
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Then, the B s — > D S ( 3 A) form factors can be expressed by 

AB.d.Ja)^ = -{m Bs -m Del{3A) ) q 3 %% Vl BM3A) tf) = ^ , 

V 2 '(q ) = (m Bs ~m Dsl(3A) )c + (q ), V 3 ' (q ) - V (q ) = c_ (q ). (18) 

2m D sl ( 3 A) 

The analysis is similar for the 1 A meson but only the w'l A terms contribute in the transition form factors, stemming 
from the structures of meson-quark-antiquark coupling vertices shown in Eq. ([8]) . 

Before closing this section, it is worth mentioning that the identification of £> s (3040) with either 3 Pi or 1 P\ 
configuration may be inappropriate. As discussed above, in the heavy quark limit the QCD Lagrangian is invariant 
under the heavy flavor and spin rotation. Thus, the heavy quark will decouple from the remanent part. One direct 
consequence is that the heavy mesons are not labeled by the quantum number 2S+1 L,j but by the total angular 
momentum si of the light quark degrees of freedom. Finite heavy quark mass corrections will break this symmetry, 
and mesons with the same spin-parity J p will mix with each other. The relation between these two bases is 

IA 1 ' 2 } = yfl'fl) - Vfl 3p i>' (19) 

As follows, we will compute the B s — > D s (3040) form factors and make predictions for physical observables in the 
relevant B s decays under these two schemes. 

III. NUMERICAL RESULTS 

Expressions for form factors shown in Eq. (|14j) involve the quark masses, hadron masses and the light front wave 
functions (LFWFs) tp' and tp' p . The latter should be obtained by solving the relativistic Schrodinger equation. However, 
in practice a phenomenological wave function to describe the hadronic structure is preferred. In this work, we will 
use the simple Gaussian-type wave function which has been extensively adopted in the literature (28l |29| 

i \ J^\ m ( P?+P'l\ dp' z e[e 2 

* = * (X2 ^ } = 4 \w) v exp (-V^J ' ^ - (20) 

Similarly for the D s (3040), as one of the 2P states, the wave function is given as 

The constituent mass for a heavy quark is usually believed close to the current quark mass. Running quark masses 
for a charm and bottom quark in the MS renormalization scheme are constrained as [39[ 

m b = (4.19t°;J*)GeV, m c = (1.27t°;{£)GeV, (22) 

while the corresponding pole masses (with one-loop anomalous dimension) are 

m b = (4.65±glg?)GeV, m c = (1.56t°;^)GcV , (23) 

with a s (niz = 91.19GeV) = 0.120 [39[ in the renormalization group evolution equation. The quark mass extracted 
from the T(1S) is 

m b = (4.67i8;5§)GeV. (24) 

In a phenomenological approach such as the quark model used in this work, the constituent quark mass is chosen 
close to the above results. 

As well known, the constituent mass of a light quark largely deviates from its current mass due to non-ncgligiblc 
contributions from gluon degrees of freedom. In this case, only phenomenological values can be employed, and in 
principle, they can be constrained within a global fit under a specific framework. In this work, we will employ the 
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values which are consistent with the previous studies in the covariant LFQM. In particular, the same values as in 
Ref. [28[ will be adopted 

m c = (1.4 ± 0.1)GeV, m b = (4.64 ± 0.2)GcV, m s = (0.37 ± 0.05)GeV, (25) 

where a comprehensive study on D, B transition form factors has been performed and many predictions are found 
in accordance with the experimental data. It is worth mentioning that Ref. (28j only gives the central value for the 
quark masses. But here we also include sizable uncertainties to estimate the sensitivity. The above values in Eq. (|25|) 
are also consistent with the results derived from the variational principle in the conventional LFQM poj . 

The shape parameter /3', which describes the momentum distribution, is fixed by the meson's decay constant, of 
which the analytic expression is derived as 



fv = TTTTZa / dx * d P± 



miTO 2 — m\ + X\Mq — p'l H —(m 2 + mi) 

Wy 



(26) 



M'4tt 3 J " ^ ± x 1 x 2 (M /2 - Mf) 

For the B s meson decay constant, we use the latest Lattice QCD result |4lj : 

J Bs = (231 ± 15)MeV. (27) 

The D s decay constant can be measured by the semileptonic D s — > pa)^ decays, for which a number of measurements 
are available. The averaged value by the HFAG at la level will be used in this work [42| 

f Ds = (256.9±6.8)McV. (28) 

At present, no direct measurement of the D* decay constant is available. Fortunately, the heavy quark symmetry 
(HQS) implies that this quantity can be related to fu s , i.e. 



f D * = J—fD M = (248.0 ±6.6)MeV. (29) 
V m D * 

Symmetry breaking effects emerge as power corrections in Aqcd I^Q (here Aqcd is the hadronic scale while tuq can 
be chosen as m c or mD s ,D*), of which the value can amount to 20%. In the present work, we will use the value 
derived from the HQS relation and neglect the power corrections. Using the above decay constants, we fix the shape 
parameters in the LFWFs as 

f3 Bs = (0.6224 ±0.0339)GeV, p Dt = (0.421 ± 0.007)GeV. (30) 

With the instantaneous Bethe-Salpeter method, the decay constant for radially excited 1 + states has been studied 
by Ref. [H, i-e. 

fo^Fi) = 204MeV, / fls(2 i Fl) = 50MeV, (31) 

but no uncertainties can be found in Ref. [l6|. With the value for /j3 s (2 3 p 1 ) and identifying D s .(3040) as one 2 3 Pi 
state, we find that the shape parameter is close to 0.3 GeV. As a first step to proceed, we will employ this value with 
an uncertainty as follows, 

/3/M3040) = (0.300 ± 0.015)GeV. (32) 

This choice corresponds to the decay constant 

/d.(2»a) = (198^°)MeV, f De ^ Pl ) = (58±5)MeV. (33) 

For later convenience, the decay constants of D~ extracted from the experimental data of D~ — > fJ.~v, and D*~ 
from HQS are also listed: 



f D = (216.6 ± 17.2)MeV, f D . = J— -f D = (209.1 ± 16.8)MeV. (34) 

771 ]j * 
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FIG. 2: The g 2 -dependence of the B s — ► D* a form factors 



Decay constants of the light mesons are obtained by making use of the experimental data of P — > Iv and r — > V v T 
decays [39| 

f„ = O.I31GeV, f K = 0.16GeV, f p = 0.2212GeV, f K * = 0.2068GeV. (35) 

It is worth pointing out that because of the condition q + = imposed in the course of calculation, form factors 
are known only at space-like momentum transfer q 2 = —q\ < 0. On the other hand, only time-like form factors 
are relevant in physical decay processes. Namely, in the exclusive non-leptonic decays, the form factor at maximally 
recoiling (q 2 ~ 0) is required; while the g 2 -dependent behavior in the full q 2 > region is needed in scmilcptonic B s 
decays. In order to have the behavior in the whole q 2 region, we adopt the parametrization 

F{ - q2) = l-aq 2 /m 2 Ba ( '+b(q 2 /m 2 Bs ) 2, (36) 

where F denotes any generic form factors V, ^0,^1,^2 for B s D*, or A, Vo, V\, V2 for B s — > D s i. These parameters 
(F,a,b) will be fitted in the space-like region (— 20GcV 2 < q 2 < 0). Then, the form factors will be extrapolated to 
the time-like region. 

With the above input parameters, our results for the form factors are collected in table HI where the uncertainties 
are from the shape parameters of the LFWFs and the constituent quark masses. Our strategy for the uncertainties is 
as follows. With the central value for quark masses, we first fix shape parameters so that these values can reproduce 
the adopted decay constants. Secondly, we will also compute the errors caused by quark masses but keeping the 
shape parameters unchanged. The final uncertainties for the form factors and physical obscrvablcs will be added in 
quadrature. We plot the g 2 -dependence of the form factors in Fig. [2}^ [6] for B s -> D*, D s i( 3 Pi), Dgi^Pi), D sl {Pl' 2 ), 

3 /2 

and D s i(P 1 ), respectively. 

Several remarks are given in order. Firstly, the B s — > D* form factors are related to the B D* form factors by 
the flavor SU(3) symmetry. Compared with the computation in Ref. [28l ] 



l-1.29g 2 /m| + 0.45(? 2 /m|) 2 ' yH ' 1 - I.30g 2 /m| + 0.31( g 2 /m|) 2 ' 

^63 A BD '(^ = — 

l-0mq 2 /m 2 B + Q.02(q 2 /m 2 B ) 2, 2 W > 1 - f .Uq 2 /m 2 B + 0.52( g 2 /m|) 2 ' 
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FIG. 3: The (^-dependence of the B s -¥ D s (3040) form factors with 2S+1 Lj = 3 Pi 




FIG. 4: The g 2 -dependence of the B s -¥ D s (3040) form factors with 2S+1 Lj = 1 Pi 



we can see that our predictions are compatible with their results. The small differences originate from the SU(3) 
symmetry breaking effects in the constituent quark, hadron masses, and decay constants of the initial and final 
mesons. 

From the results in table HI we can see that although most B s — > D s i form factors are smaller than the relevant 
B s — > D* ones, they are still large enough. It implies that the B s — > D s i decay channels are very likely to be 
observable. In particular, the large form factors yB» D si( p i ) anc j yB,D al ( p ± ) ^ TCC ^y predict the sizable production 
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< -0.2 




FIG. 5: The q -dependence of the B a — > D s (3040) form factors with the quantum numbers P± 




FIG. 6: The q -dependence of the B a — >• D s (3040) form factors with the quantum numbers P[ 



3/2 



rates in nonlcptonic B s decays, which will be shown in the following section. 

One important feature of the B s — > D* form factors is that they become larger with the increasing q 2 . However, 

not all the B s — > D s \ form factors have this universal behavior. For instance, the form factor yB» D »n p i) decreases 
with the increasing q 2 . 

Most form factors have the uncertainties of the order 10% but some of them can reach 20%. The uncertainties are 



denoted by the bands in Figs. (3-6). The q -dependent parameters of V 2 



B S D S {A 3 ' 2 ) 



also suffer from large uncertainties. 
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TABLE I: B s -t D* and B a 
of quantum numbers. 



Z) B (3040) form factors. Four sets of results are shown for D s (3040), with different assignments 



V BsD* 
A 

A B °°> 



DM 3 Pi 
i( 3 Pi 

D s i( 3 Pi 
1(^1 
li'Pl 



A 
V 

yB. D 



j4 B 3 i3 sl (F 1 



T7T 



V, 



B„D sl (Pl> 



B s D al (Pl< 



yB B D sl (P^ 
A B s D sl (p3> 



V, 



I) 



sl(Pf 



B s D sl (P*' 



sl(Pf 



0.74 
0.63 
0.61 
0.59 



+0.05 
0.05 
+0.04 
0.04 
+0.03 
-0.03 
+0.04 
0.04 



+0.03 
0.03 
+0.02 
0.02 



07+0-09 
— 0.10 



0.29 
0.27 



0.07 
0.12 
0.52 
0.62 
-0.20 



+0.02 
0.02 
+0.01 
0.01 
+0.06 
0.06 
+0.06 
0.06 
f0.03 
-0.04 



n , 7 +0.03 
u - i '-0.02 



0.08 



+0.05 

0.05 

+ 0.11 



-0.43 _ 

n i o+0.03 

n o7 +0.02 
u ' z ' -0.02 

+0.04 



0-57; 

106 +0.04 

-0.09 



0.04 
0.06 



1.27 
1.05 
0.82 
0.97 



+0.13 
0.12 
+0.06 
0.06 
+0.05 
0.05 
+0.09 
0.08 



+0.05 
0.05 
+0.02 
0.02 

0.90_g;H 

+0.02 
0.02 
+0.02 
0.02 
74+0-10 



0.39 
0.37 



0.07 
0.16 



0.68 
-0.29 



+0.07 
0.07 
+0.06 
0.06 



-0.22 



+ 0.04 
0.04 

O-^toTo 
-0.37_°^ 
_n oq+o.05 

u - zo -0.05 

S6+ 03 

U.OO_o.o3 

+0.07 
0.07 

+0.05 
0.06 

+0.03 
0.04 



0.82 
1.06 
-0.03 



(i, 



1.34 
1.30 
0.71 
1.24 



+0.11 
0.10 
+0.09 
0.09 
+0.09 
0.08 
+0.10 
0.09 



1.36 
1.52 
-0.35 
-0.25 
1.50 
1.73 
0.46 
2.03 



+0.12 
0.12 
+0.10 
0.10 
+0.15 
0.18 
+0.39 
0.46 
+0.10 
0.10 
+0.11 
0.11 
+0.12 
0.12 
+0.12 
0.12 



1 29+ - 14 

1 - Z9 -0.15 

9 i c+0.08 
z - ld -0.25 

-o.87tg:g 

i 9 i+0.16 

^^-o.io 

1 41+0-10 
- L - 4 - L -o.io 

1 67+ - 11 
L - 01 -0.11 

01+ 07 

u - ul -0.08 
"-10 



0.57 
0.54 
0.08 




+0.15 
0.13 
+0.14 
0.12 
+0.05 
0.04 
+0.13 
-0.11 



0.26 
0.30 
0.33 
0.53 
0.71 


0.12 
1.94 



+0.06 
0.04 
+0.06 
0.05 
+0.14 
0.11 
+0.31 
0.23 
+0.17 
0.14 
+0.09 
-0.08 
+0.04 
0.03 
+0.32 
0.26 



0.14 
1 



+0.06 
0.03 

+3.7 
-1.0 



R7+0- 38 
U - D ' -0.26 

21+ 07 

U- zl _0.04 

40.+ 008 



0.42Z 
0.17 

2o 



+ 0.05 
0.04 

+ 17 
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It is fortunate that the magnitude of this form factor is small. Thus, it has little influence on physical processes. 



IV. APPLICATIONS TO B s DECAYS 



In this section, we will use the above determined form factors to predict semilepotnic B s — > D*(D s (3040))£p and 
nonlcptonic B s decays. For nonlcptonic B s decays, the factorization assumption will be used to decompose the decay 
amplitudes into the form factors and the decay constant of the emitted meson. In particular, we will only consider 
the color-allowed decay channels where factorization assumptions work well and the nonfactorizable contributions are 
typically small. 



A. Semileptonic B s decays 

With the form factors available, we can investigate the semileptonic B s —> D* s and B s —> D s (3040) decays 

dv M v drf_ 

dq 2 ~ V dq 2 ' 

i— L,± 

dr% ± \G F V cb \ 2 ^ ( m 2 \ 2 M M 

~W ~ 256m%yq 2 \}~~?) (Xi ' X±) (38) 

where M denotes D* s or D s (3040), A M = A(m| s , m 2 ^, q 2 ), and X(a 2 ,b 2 ,c 2 ) = (a 2 - b 2 - c 2 ) 2 - 46 2 c 2 . The subscript 
(L, ±) denotes the three polarizations of the vector and axial-vector cs meson along its momentum direction (0, ±1). 
In terms of the angular distributions, we can study the forward-backward asymmetries (FBAs) of lepton which are 
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defined as 



dA% Sldz{dT M /dq>dz)- j\dz(dT M /dq 2 dz) 



dq 2 



Jo dz(dT M /dq 2 dz) + dz{dT M /dq 2 dz) 



where z = cos 6 and the angle is the polar angle of the lepton with respect to the moving direction of D*(D s i) in 
the lepton pair rest frame. Explicitly, we have 



dA 



FB 



1 



dq 2 



dA 



FB 



dq 2 



x D L ° + X? + x D J 



x D sl + x D 3l + X D S 



2m 2 n /X c Zh^(q 2 )A (q 2 ) - ^.fh^A^Viq 2 ) 



where 



d: 



X 



X? = 



2 r 
3 

2q 2 



(2q 2 + m 2 )(h°° (q 2 )) 2 + 3X D *m 2 A 2 (q 2 ) 



(2q 2 +m 2 ) 



(m Bs + m D ,)Ai(q 2 ) qp 



m Bs + m D , 



-V(q 2 ) 



X 



1 



2 



{m 2 B - m 2 D , - q 2 )(m Bs + m D *)A 1 (q 2 ) - 



(2q 2 + mj)(h^ (q 2 )) 2 + 3X Dsl mjV 2 (q 2 ) 



2 -^(2q 2 + m 2 ) 



m Bs + m D , 



-I 2 



-A 2 (q 2 



(m Bs - m Dsl )Vi(q 2 ) T A(q 2 ) 



1 



2m Dsl 



(m 



3l - <l 2 )( m B s - m Dsl )Vi(q 2 ) 



X 



D , 



m Bs - m Dsl 



-V 2 (q 2 ) 



r A/ = rf + rf + r M , A% 



i 



dq' 



sign(z)dz(dT J dq dz) 



with 



(iriB s —mm) 



dq 



rlT M 



dq 2 



Since there are three different polarizations, it is also meaningful to define the polarization fraction 



r L + r+ + r_ 



(39) 



(40) 



Integrating over the q 2 , we obtain the partial decay width and integrated angular asymmetry for the B s — > D s j decay 
modes: 



(41) 



(42) 



Our results for the branching fractions, FBAs and polarizations are given in Table HJ where V c b = 0.0412 is 
employed (39j . The uncertainties are from the form factors. Several remarks are given in order. Our result for 
the B(B S -> D*W) is well consistent with the QCD sum rules in Ref. B(B S -> D*W) = (5.6 ± 0.2)%, but is 
smaller than the prediction in Ref. [H by roughly 40%, B{B S -> D*ei>) = (7.09 ± 0.88)%. At present there are no 
experimental data on the semileptonic B s decays into D*. Interestingly, our predictions are close to the B — > D*lu 
data [39[ which can be related to the B s decays by the flavor SU(3) symmetry 



B(B° -> D* + l~ v) = (5.16 ±0.11)%, B(J3°->i3* + r-p r ) = (1.6±0.5)%. 



(43) 



Since the flavor symmetry breaking is at a few percent level, our results are also supported by the nonleptonic B and 
B s decay data [39| 



B(B° -> D+7T- 



(2.68 ± 0.23) x 10" 3 - B{B° S -> D+Tr") = (3.2 ± 0.9) x 10" 



(44) 
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TABLE II: Branching ratios, polarization fractions and angular asymmetries for B a — > D't, and B 3 — s> D s (3040) 





e = n 


t = T 


B(B° -»• D* a +£is e ) 
f L (B°^D* + ev e ) 
A FB (B°s -»• D'+e.vt) 


(5.2«; 6 6 ) x 10- 2 
(51.311;?.) x 1(T 2 
(-21.711;*) x 10- 2 ) 


(1.3l°; 2 ) x lO" 2 
(44.5±g-|) x 10~ 2 

(-s.eil;!) x 10- 2 


B(B° a -+ D+^PJlut) 
f L (B° s D+fP^ivt) 

A FB (B°^ D+fPjm) 
B(B° S ^D+CPi)£u e ) 
f L (B° s ^D+CPi)^i) 

A FB (B° S ^ DjXPi)^) 


(2.3±g-|) x 10 -3 
(43.3±t"|) x 10- 2 
(-39.7±*;g) x lO" 2 

(2.9l°;«) x 1CT 3 
(84.01?,; 2 ,) x 10~ 2 
(-6.811;*) x lO" 2 


(6.1±1;|) x 10" 5 
(36.0t£;lj) x KT 2 
(-16.6itS) x 10" 2 

(5.2ll;g) x 10" 5 
(65.91 3 ;;) x 10- 2 
(19.0l 2 ; 5 6 ) x lO" 2 


e(B s -> fji(Pi )ivt) 

A FB (B° -► D+(Pl /2 )£De) 
B(B° a -> D+(P? /2 )£v t ) 
f L (B° s ^ D+ipVyvt) 

A FB (B° -> D+(P? /2 )£u t ) 


/ o r + 1 . 1 \ m-4 

(3.5l 10 ) x 10 
(9.81* 1 /) x 10" 2 
(-65.0lJ 3 2 4 ) x lO" 2 
(4.0lg; 4 x 10" 3 
(63.9llo) x 10~ 2 
(-22.2I2;?) x 10 -2 


(9-9-3.5) x 10 
(18.0l 2 -j) x 10~ 2 
(-47.6l^|) x 10" 2 

(9.rt8;l) x 10- 5 

(49.911-g) x 10~ 2 

(2.0l° o ;D x 10- 2 



Furthermore our results show that the polarization fractions and angular asymmetries are sizable and can be tested 
at the hadron collider and the Super B experiment in the near future. 

Compared with B s — > D*W, the BRs of the B s — > D s (3040)lP (/ = fx) decays are typically smaller by one order 

1/2 

of magnitude. However, if D s (3040) has the quantum number P 1 , the BR would be suppressed by two orders of 
magnitude as a consequence of the tiny form factors. The mass of the D s (3040) is larger than that of D* , making 
the phase space in B a — > D s \tD t much smaller and inducing further suppressions of the BRs. Nevertheless, the BRs 
are still large enough to measure in the future. For example, the LHCb experiment will produce more than 10 6 B s 
mesons per running year, implying that a plenty of the _D S (3040) events will be generated. 

Although the production rates for Z} s (3040) under the two different descriptions 3 Pi and l P\ are similar, their 
polarizations are quite different. For the 3 Pi assignment, the form factor V\ is much larger than the other ones. The 
interference between different form factors does not manifest. Since both longitudinal and transverse polarizations 
depend on this large V\, the polarization fraction is close to 50%. For 1 P\, the magnitude of V2 is enhanced and 
it has a different sign with V\. These two form factors will provide constructive contributions to the longitudinal 
polarization as shown in Eq. (j40|) . and lead to a larger value for fjj. 

The angular asymmetries in B s — > D s \lv with D s \ = ( 3 Pi, 1 Pi) are also different. If the lepton is a light muon, 
the mass terms in the angular asymmetries as in Eq. ([59]) are negligible. Thus, the form factor product V\(q 2 )A(q 2 ) 
determines the size of the angular asymmetries. For 3 Pi, the product is much larger than that of the 1 P\ assignment. 
So is the angular asymmetry, but with a negative sign. If the lepton is changed to a massive tau, the two terms in 
Eq. (|39[) give destructive contributions, which will reduce the modulo of the angular asymmetries for 3 Pi and change 
the sign for the 1 Pi. 

The above analysis is also similar for the other assignment basis. One exceptional case is the polarization fraction for 

1 /2 

the Pj . All form factors are relatively small, and the cancelation between V\ and V2 in the longitudinal polarization 
results in a small polarization fraction 



B. Nonleptonic B a decays 

After the contraction of the highly-virtual degrees of freedom such as the W, Z bosons, the effective Hamiltonian 
governing nonleptonic B s — > D s ,j decays is given as 

H eS = ^j=Y. V - bV M c i i + c 20 2 ) + h.c., (45) 

p,Z> 

with the four quark operators 

Or = [c° 7 M (l - J 5 )b a ][D^(l ~ l5 )p% 2 = - 7»)^][^7m(1 " 7sK], (46) 
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and C\ and C2 being their Wilson coefficients which arise from the perturbativc coefficients over the rrib scale. We 
use ai = C\ + C2/3 = 1.07 [34j]. In the above equation, it is labeled p = u,c and D = d,s, while a and /3 are color 
indices. The CKM matrix elements are employed as V ud = 0.97418, \V cd \ = 0.23, V cs = 0.997 and V us = 0.2255 j^. 

Nonleptonic B s — > D s jM decays involve one more meson and the two final state mesons will get entangled with 
each other. In case of the M being a light meson, it is feasible to show that the decay amplitude can be decomposed 
into two individual parts. This method is known as factorization [4614481 ] . The light meson moves very fast in the B s 
rest frame and it is made of collinear quark fields and gluon fields. On the contrary, the B s meson and the recoiled 
D s j arc almost at rest and the dynamic degree of freedom is the heavy quark, the light antiquark with quantum 
fluctuations. The complete set of the degrees of freedom in these heavy mesons has typically small momentum 1 and 
is characterized by the soft fields. Interactions between the two different sectors are forbidden at the leading power 
and these two sectors are dynamically separated. The elegant proof rooted in the quantum field theory has been given 
in Ref. [49|, where an effective field theory built in Ref. [50[ has been used. 

Under such an effective theory, the matrix element of the four-quark operators in two-body B s decays is expressed 
in terms of the B s — > D s j form factors, which have been computed in the above, together with the convolution of a 
hard kernel T(u) and the light-cone distribution amplitude <p(u) of the light meson. The hard kernel can be computed 
in the expansion of a s , and particularly at the leading order T(u) = 1 + 0{a s ). In this case, the analysis based on the 



effective theory recovers the naive factorization in which the decay amplitudes for B & 
are given as 



DZ 



and B, 



D 



*+, 



A(B° S -+ D*+n- 



A L (B° S -+ Dl+p-) 



A N (B° S -> D*+p- 



A T (B° S -> D*+p- 



G f I 

-JLVcbV^axU A (ml)^J\(m? Bs , m%. , m|), 

—j=-V cb V ud aif p h (m p ), 
-iG F 

-%G F 



V c bV* d aif p (m Bs +m D *)m p A 1 (m 2 p ), 



VcbV* d a!fp 



\(m Bs ,m 2 p ,m 2 D *) 



-m p V(m 2 ). 



V 2 \ m B s + m D , ) 

In the above, the amplitude of B s — > D* + (Pd*) P~ (P p ) has been decomposed according to the Lorentz structures 



(47) 



A = A L + (* Dt (T) ■ e* p {T)A N + iA T t a p 1P ^f 



2Pl,PP p 



X(m Bs ,m 2 D ,,m 2 p ) 



(48) 



We will also consider the decay channels in which the light meson (it , p ) is replaced by A" , D ,D S 
K*~ , D*~ , D*~. The partial decay width of B s — s- D*P, where P denotes a pseudoscalar meson, is given by 



T(B S -> D*P) 



\P\ 



\a(b s ^d;p)\ 



(49) 



with \p\ being the three-momentum of the D* in the B s meson rest frame. For B s 
is the summation of the three polarizations 



T{B S ^D* S V) = 



8irm B 



D*V , the partial decay width 



Ao(B s ^D*V)\ +2\A N (B,-^D*V)\ +2\At(Bs^D*V)\ 



(50) 



In the decay amplitudes for the channels involving D s i, the form factors are replaced correspondingly and ms s +mn* 
is replaced by ms 3 — fnD 3l - 

Before proceeding, it is worth mentioning that the proof of factorization of B s decays into two charmed mesons 
such as B s — > D*D* is absent in the literature. Nevertheless, since only the color-allowed tree diagrams are considered 
in this work, the nonfactorizable diagrams are typically small and the factorization usually works very well (46l - |48| . 
For instance, our predictions of the branching fractions of these channels are well consistent with the experimental 
results. This verifies the applicability of the factorization approach. 



For the heavy quark, it means the residual momentum after separating the large label momentum. 
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TABLE III: Branching ratios of nonleptonic B a decays into D* s and a comparison with theoretical results [5 
mental data 



551 ] and experi- 



B° - 


-+ D* a 


f 7T~ 


2.8 x 10 


m- 






2.1 x 10 


m- 


-> i>r 




3.3 x 10 


m- 






7.4 x 10 


m- 




V 


8.9 x 10 


m- 






4.8 x 10 


m~ 




D*~ 


1.0 x 10 


m- 




or 


2.9 x 10 



Factorization [53] 



-3 
-4 



perturbative QCD [54. 551 



(2.42 
(1.65 
(2.7 
(7.0 
(5.69 
(3.47 
(3 

(9.9 



+ 1.37 
1.06 
+ 1.06 
0.82 
+ l-9\ 
1.4) 
+4.8\ 
3.8) 
+3.59 
2.84 
+2.24 
1.72 
+2.9\ 
-2.3) 
+7.7 \ 
6.2) 



I x 10" 
i x 10" 
x 10"' 
x 10 _: 
i x 10" 
i x 10" 
x 10"' 
x 10 _: 



This work 



(3.5t°' 4 ) x 10- 3 
(2.8t° 3 3 ) x lO" 4 
(3.7l°;l) x 10- 4 
(9.2± 4 ; 4 ) x 10~ 3 
(ll.Stf.;?) x IO" 3 
(5.5t°;») x 10- 4 
(8.6+ 4 ;° 9 ) x lO" 4 
(23.6l 4 ;° 8 ) x 10- 3 



Exp. 



(2-4±S:S) x io- 



(ll.Sli; 4 ) x 10- 



(3-l±^)% 



TABLE IV: Branching ratios of nonleptonic B 3 decays into D a (3040) with different assignments 





TV 


K~ 


D~ 


D- 




(3.2±° B 5 ) x lO" 4 
(1.2l«; 3 ) x lO" 3 


(2.5±°;1) x IO" 5 

(9-3±2;g) x io~ 5 


(1.4±°;») x IO" 5 
(5.7t\i) x 10 -5 


(2.8«; 4 ) x IO" 4 

x io- 3 


B° a ^D+(Pl /2 ) 
B° S ^D+(P? /2 ) 


(3.0t|;|) x 10" 5 
(1.5lg;5) x lO" 3 


(2.3±|;1) x 10~ 6 
(1.2±°;?) x IO" 4 


(1.5±?;2) x IO" 6 
(6.9±i;J) x 10~ 5 


(3.0+ 4 - 4 ) x IO -5 
(1.4i™) x IO" 3 




P~ 




D*~ 


d;- 


B°^D+( 3 Pi) 
B° B ^D+CPx) 


(i-3±S;|) x io- 3 

(3.0ig;S) x IO" 3 


(6.4±°;J) x IO" 5 
(1.4l°; 3 ) x IO" 4 


(6.6±i^) x IO" 5 
(4.5ir 8 ) x IO" 5 


(1.4±g;l) x IO" 3 
(9.0t?;°) x IO" 4 


B° a ^D+(P? /2 ) 


(1.6t£;°) x IO" 4 
(3.3j£?) x 10~ 3 


(8.5l 4 - 4 ) x 10~ 6 
(1.5i° ; 3 ) x 10- 4 


(l.ltg;5) x IO" 5 
(8.9±g;£) x 10~ 5 


(2.4ti'J) x IO -4 

(i.9± :l) x io~ :i 



Using the above parameters, we directly obtain the branching fractions which are collected in Table Hill and Table HVl 
For comparison, we also collect two sets of different theoretical predictions on the B s — s- D* decavs~ |53l - [55| and the 
relevant experimental data [13, EH ■ From Table Mil we can see that our predictions on B s — s- D* s decay channels 
are consistent with the experimental data within uncertainties. This may also imply that our results for B a — > 
£> s (3040) are reliable. Results in the factorization method [53| are also close to our results but the perturbative 
QCD predictions [U, [55| are typically smaller. For instance, the central value of the perturbative QCD result of 
B° s -> D* + D*~ is smaller than our result by a factor of 2.5. 

Our uncertainties are from the form factors, i.e., from the quark masses and the shape parameters. Uncertainties 
in the perturbative QCD approach shown in Table HTT1 are from three sets of input parameters: (1) decay constants of 
f B , = (0.24 ± 0.03) GcV (in Ref.[H|) or f Bs = (0.23 ± 0.03) GeV (in Ref.fH), and the shape parameters of the B s 
wave functions uif, = (0.50 ± 0.05) GeV; (2) the factorization scale (from 0.75t to 1.25t not changing the transverse 
part) and the hadronic scale Aqcd = (0.25 ±0.05) GcV; (3) the CKM matrix elements V cs and V u d- Among them, the 
hadronic inputs are found to give the largest uncertainties, while the ones from CKM are the smallest. For instance, 
the branching fraction of B® — > D* + ir~ with these three kinds of errors is given as 

B(B° -> D:+0 = (2A2tl- 1 7 ltrr^l) x (51) 
Uncertainties in the PQCD approach are larger than our results for several reasons. 

• Compared with the Jb s adopted in this work, the uncertainty in their computations is larger by a factor of 2. 

• In the PQCD approach, the B s wave function has such a property that the decay constant can be factorized out 
and then fg, is only the normalization constant of the wave functions. In the PQCD approach, the authors have 
independently estimated these two uncertainties and added them together as their final results. On the contrary, 
in the covariant LFQM the decay constant can not be factorized, but has been expressed as a convolution form 
as shown in Eq. ([2"o]) . In this case, only the decay constant is the origin for the uncertainties, which will be 
smaller. 



• The uncertainties caused by factorization scales and hadronic scales in the PQCD approach characterize higher 
order QCD corrections and can be viewed as the model-dependent errors (or systematic errors). They are 
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smaller than but still in similar magnitude with the ones from the B s wave functions. One particular type of 
model-dependent uncertainties in this quark model concerns the zero-mode contributions. A direct evaluation 
of the form factors in the light-front quark model suffers from the non-covariance problem. It is resolved with 
the inclusion of the zero-mode contributions, i.e. the Z-graph as in the conventional LFQM, and the direct 
treatment of the spurious terms with the replacement rules as in the covariant LFQM. Different results can be 
produced. However, at present such computation in the conventional light-front quark model is not available. 
Therefore, it is not possible to estimate the systematic errors. 

In B s — > D s i decays with a pseudoscalar meson emitted, the decay amplitudes only involve the form factor Vq. The 
large form factors for the P^ 2 and l P\ result in large BRs. Similarly, the large uncertainty for Vq{B s — > D s i(P^ 2 )) 
also gives large uncertainties to the BRs. 

Channels with the emissions of K~ and D~ arc suppressed by roughly one order of magnitude compared with the 
case emitting a pion and D s meson. This is due to the hierarchy in CKM matrix elements: \V US \ ~ \V c d\ ~ 0.22 <C 
\V ud \ ~ \Vc\ ~ 1. 

Under different quantum number assignments for D s (3040), the branching fractions B s — > D s (3040)p and B s — > 
D S (304:0)D* can reach the order of 10~ 3 . But the modes with the emissions of K* and D* are smaller by 1 to 2 orders. 
Nevertheless, all of these channels have a promising prospect on the LHCb experiment with sufficient statistics for 
the B s production. 

V. CONCLUSIONS 

Using the covariant light-front quark model, we have studied the B s — > Z) s (3040) form factors, including the 
B s — > D* form factors as a byproduct. The form factors are computed in the space- like region and extrapolated 
to the physical region through a three-parameter form. Those form factors are used to predict the semileptonic 
and nonlcptonic B s decays. In particular, the predictions on branching fractions, longitudinal polarization fractions 
and the angular asymmetry, arc provided. Our results for B s — > D* decays are consistent with the experimental 
data for the B s and B decays which are related by the flavor symmetry. For B s — > _D S (3040), we find that the 
branching fractions arc large enough to be observed on the LHCb experiment. For instance, the BRs of semileptonic 
B s —> D s (30A0)lP, nonlcptonic B s — s- D s ip, and B s — s- D s iD* can reach the order of 10~ 3 . In spite of similar 
production rates, the different assignments for the D s (3040) can be distinguished by the polarization fractions and 
angular asymmetries. We expect more experimental results from the B factories and LHC etc in the near future. 
They would be helpful to provide deeper insights into the cs spectroscopy. 
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Appendix A: Some specific rules under the p~ integration 

When performing the p~ integration one needs to includes the zero-mode contribution to solve the noncovariance 
problem. This amounts to performing the integration in a proper way in the covariant LFQM (27l. [28|. In particular 
for p' lfl and p'iuPi v we have 

&,&u = <W4 2) + P»PvA (2) + {P^q v + q^P v )A {2) + q^q v Af\ (Al) 
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with the symbol = denoting that these equations are valid after integration. are functions of Xi s 2, p' 2 , v'\_ ' <7-L 
and q 2 



Z 2 = N[ + mf - ml + (1 - 2x{)M' 2 + {q 2 + q ■ P)- 



<l 2 



4(1) _ x i A W _ 4 (i) g± • ^ 4(2) _ _* (p!l • 



9" 9- 
4 2) = A^A^\Af = {A^f-^Af\ (A2) 
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